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The cross sections for a variety of diffractive processes in proton-nucleus scattering, associated with 
large gaps in rapidity, are calculated within an improved Glauber-Gribov theory, where the inelastic 
shadowing corrections are summed to all orders by employing the dipole representation. The effects 
of nucleon correlations, leading to a modification of the nuclear thickness function, are also taken 
into account. Numerical calculations are performed for the energies of the HERA-B experiment, and 
the RHIC and LHC colliders, and for several nuclei. It is found that whereas the Gribov corrections 
generally make nuclear matter more transparent, nucleon correlations act in the opposite direction 
and have important effects in various diffractive processes. 
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I. INTRODUCTION 

In hadron-nucleus collisions at high energies nuclei act 
almost like "black" absorbers. Therefore the optical anal- 
ogy should be relevant and diffraction appears to be an 
important process. Experimentally diffraction appears 
as large rapidity gap events, when the debris of the pro- 
jectile hadron and the nucleus occupy only small rapid- 
ity intervals close to the rapidities of the colliding par- 
ticles. The optical analogy is employed by the Glauber- 
theory [l[ of hadron-nucleus interactions, which assumes 
additivity of the scattering phases on different bound nu- 
cleons. This is a single channel approximation assum- 
ing that absorption, i.e. inelastic interactions, generates 
via the unitarity relation only elastic scattering. In real- 
ity, diffractive excitations of hadrons frequently happen, 
and the Glauber approach was generalized to a multi- 
channel case by Gribov 0. The corresponding correc- 
tions to the Glauber approximation are known as inelas- 
tic shadowing, or Gribov corrections. Unfortunately, the 
multi-channel problem needs detailed experimental in- 
formation, which is mostly unknown. One has to know 
all diffractive amplitudes, diagonal and off-diagonal, for 
different diffractive excitations of the hadron. Even the 
lowest order correction contains an unknown attenuation 
factor for an excited state propagating through the nu- 
cleus m. 

One can sum up the Gribov corrections to all orders by 
switching to the interaction eigenstates which were 



identified in Q as color dipoles, and where the dipole 
approach to high energy collisions was proposed. This 
phenomenology needs lesser input and the key ingredi- 
ent, the dipole- nucleon cross section, is flavor indepen- 
dent and can be studied in different processes. 

This method can be applied also to lepton- or photon- 
nucleus collisions [1, 0,11], where leptons and photons 
display hadronic properties. A detailed study of the in- 
elastic shadowing corrections to different diffractive chan- 
nels in proton-nucleus collisions was performed, within 
the dipole approach, in 0, • 

Here we are going to enhance the accuracy of the calcu- 
lations presented in 10], by improving the model for the 
nuclear wave function. Namely, most of calculations for 
nuclear shadowing effects have relied so far on a simplified 
model of an uncorrelated single particle density distribu- 
tion in the nucleus. This model in particular ignores the 
well known experimental evidences for the existence of a 
strong repulsion core between nucleons. Such a repulsion 
should lead to short-range NN correlations in the nu- 
clear density function, which in turn should modify the 
effective nuclear thickness function controlling diffractive 
processes. 

The consideration of possible effects from nucleon- 
nucleon (NN) short range correlations (SRC) appears to 
be particularly interesting, in view of recent experimen- 
tal data on lepton and hadron scattering off nuclei at 
medium energy, which provided quantitative evidence on 
SRC and their possible effects on dense hadronic matter 
[ll| . Moreover, a recent calculation of the total neutron- 



nucleus cross section at Fermilab energies has indeed The main problem in evaluating nuclear effects is there- 
shown relevant effects from SRC even at high energies fore the choice of the nuclear wave function. ipo(l, . . . , A). 

m. 



II. GLAUBER FORMALISM 

The key assumption of the Glauber model is that the 
hadron-nucleus partial elastic amplitude at impact pa- 
rameter b has the eikonal form 



rt> A (b;{T j , Zj }) = i-l[ [i-r*" (b-T k ) 



(i) 



fc=i 



where {lj,Zj} denote the coordinates of an i-th target 
nucleon; iT pN is the elastic scattering amplitude on a 
nucleoli normalized as 



of" =2 J d 2 bReT pN (b) 



(2) 



Further, one should calculate the matrix element of 
the amplitude |T]) with the nuclear wave function, 
4>o(r\, f2, . . . , r A ) — ipo({rj}) = |0). Here we introduce 
new notations, 



G A (b) = (0^(6; {^})|0) = 1- (0|JJ G pN (b-h)\0) 

i=i 

(3) 



where 



G pN (b-h) = l-T pN (b-h) 



A. Single particle approximation for the nuclear 
wave function 



The most popular model for the square of the nuclear 
wave function appearing in the Glauber formalism is the 
approximation of single particle nuclear density [3t| . 



^(fi^.^DI^JIpit^zi) (5) 

3=1 



where 



PA 



A 

(h,zi)= I Y[d 3 n\y A ({rj})\ 



(6) 



Within such an approximation the matrix element be- 
tween the nuclear ground states reads, 



1 



I --- / d 2 lT pN (l) I dzp A (b-l,z) 



•(7) 



(4) Correspondingly, the total pA cross section has the form, 



Jd 2 b J 1 - 1 - j J d 2 lT pN (l)T A (b- f) 



>!:.) = 2 Re / ,l 2 h i I - 

2 / irh { I cx|. 



>,Jd 2 bx |] 



■\a^ t {l-ia pN )Tl{b) 



r 



where a P N is the ratio of the real to imaginary parts of 
the forward pN elastic amplitude; 



pN 



and 



T A {b) 



d 2 lReT pN (l)T A (b-f); (9) 



dz PA (b,z), (10) 



is the nuclear thickness function. We use the Gaussian 
form of T pN (l), 



ReT pN (l) = 



pN 



-I 2 



(11) 



Notice that in Eq. ([5]) and in what follows we use 
the exponential approximation of large A only to sim- 
plify and clarify the formulas. For numerical calculations 
throughout the paper we always rely on the exact expres- 
sions, such as the first part of Eq. ([5]). 

The Glauber approach is a single channel model, there- 
fore it is unable to consider diffractive excitation of the 
proton. However, a part of diffractive excitation of the 
nucleus occurs without excitation of the bound nuclc- 
ons, when the nucleus just breaks up into free nucleons 
and nuclear fragments. Such events, pA — * pF, are called 
quasielastic and can be calculated within the Glauber ap- 
proximation. Summing up the final states of the nucleus 
l-F), applying the condition of completeness, and extract- 



2 



ing the contribution of the ground state of the nucleus, one gets, 
I 



pA 



F 

= [ d 2 b \{o\r pA (b)\F) f (F\r pA {b)\o) -\{o\r pA (b)\o)\ 



d z b 



T pA {b)\ 



|(o|r pA (6)|o) 



(12) 



In the first order in nuclear density the first term in I) [r pJV (7)] 2 = T^(b)a p ^ . Both terms together result in 



this expression, ( 



r pA (6)r 0), contains, besides the the exponent a? T A (b). Then the quasielastic cross sec- 



usual linear term Eq. @, the quadratic term J d 2 sT A {b~ 



tion gets the form, 



a pA ( P A P A*) = J d 2 b {exp [-<f T» 



exp 



(&)] } • 



(13) 



B. Nucleon correlations 



Equation ([5]) represents only the lowest order term of 
the square of the full nuclear wave function \ipo\ 2 . As a 



matter of fact, the latter can be written as an expansion 
in terms of density matrices [3, [l3| as follows: 



(14) 



3 = 1 



in which the single particle density pi{ri) is 

f 2 A 

Pi(n) = / \ipo(ri,f 2) ...,f A ) \ 2 ~[[d 3 r, 

J i=2 

and the two-body contraction A is 

A{n,rj) = p2{n,fj) - pi(fi) pi(fj) 

The two-body density matrix 



A 

p2(ri,r 2 )= I \ip (ri,r 2 ,... 1 f A )\ 2 Y{d 3 - 

satisfies the sequential condition, 

d 3 r 3 P2(n,fj) = pi(n) 



i=3 



(i<j)^(k<l) 



(15) 



(16) 



(17) 



(18) 



r 



which leads to the basic property of the two-body con- 
traction 



d 3 rj A(n, fj) =o. 



(19) 



Notice that the single particle density appearing in Eq. 
(TTJ} is normalized to one, so that the densities de- 
fined by Eq. © and Eq. ([15]) are simply related by 
p A {r) — Api(r). It should be stressed that in Eq. (fT4"|) 
only unlinked contractions have to be considered, and 
that the higher order terms, not explicitly displayed, in- 
clude unlinked products of 3, 4, etc two-body contrac- 
tions, representing contributions to two-nucleon correla- 
tions, and unlinked products of three-body, four-body, 
etc, contractions, describing three-nucleon, four-nucleon, 
etc correlations. We will give now a short derivation of 
the total cross section including two-nucleon correlations 
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(more details will be given elsewhere [14]), Taking into 
account all terms of the expansion ([M]) containing all 
possible numbers of unlinked two-body contractions, Eq. 



([5]) can be written in the following form which yields the 
usual Glauber profile when A = [H,[lJ]: 



.A A 
G A {b)= / l[d 3 r k \i>(l...A)\ 2 l[Gi> N (b-k) = 

J fe=l 1 = 1 

/A n A A 

Y[d 3 r kPl (r k )G pN (b-l k ) + Y, / li^rkAin,^) [] Pl {n)G pN {b - l k ) 

k=l i<j J k=l l=£i,j 

r A A 

Hd 3 r k A(n,r 3 )A(r p ,n) J] Pl {r m )CF N {b - l k ) + ■ ■ ■ 

- A _L ^ ( ** 7, 1 LI A I 



i<j^P<l k=l 



G^(b) + G^(b) + G^(b) + . 



(20) 



r 



where the superscript denotes the number of two-body where 
contractions in the given term, each term containing 
Glauber profiles to all orders. For each nucleus we have 
considered all terms of the series (|2Tfl) ; the first term, cor- 
responding to the single particle approximation of Eq. 
© being 



X(b) = / d 3 nd 3 r 2 A(fi, 7^)^(6 - h)TP N (b- f 2 )(23) 



r A 

Gf{b) = Y[d 3 r kPl (r k )GP N (b-l k 

^ k=l 



and 



I- j I d^npAin^ib-h) 



(21) 



Y(b) 



A 



(24) 



and the n-th terms 



G%\b) 



resulting from the basic properties of the two-body con- 



2 n n\(A- 2n)\ 



X n (b)Y A ~ 2n (b) (22) traction: f d 3 n,jA(fi,rj) = and 



d 3 r 1 d 3 r 2 A{r ll r 2 )G vN {b- h)G pN (b - f 2 ) = / d 3 r 1 d 3 r 2 A(f 1 , f 2 )TP N (b - \b - f 2 ). 



(25) 



Eq. (|2"0"1) can now be written as follows 



n=0 



2«n!(A-2n)! a»i [r WJ 2^ 2 " nl 1 WJ 



n=0 



(26) 



Using, for ease of presentation, the optical limit approx- 
imation 



4 



A 



(27) 



the insertion of Eq. (|2T|) into Eq. (|2^|) leads to the final 
result: 



Ga (b) — 1 — exp 
= 1 — exp 



- j &n pA(n)T{b-h)- 



- I d* ri d 6 r2 A A (ri,r 2 )T(b - h) T(b - l 2 , 



(28) 



where 



A A (n , r 2 ) = pf (n ,r 2 )~ PA {n) p A (f 2 ) . (29) 

which obviously differs by Eq. (TI75)) simply by a factor ^4 2 , 
and 



with rjf(6) given by Eq. fl9|) and 



(30) 



AT» = 



(l-ia pN ) f 2 



d^Zid^ 2 Ai(/ x ,/ 2 ) 



where 



x Rer pJV (fo-fi) Rer pJV (6-r 2 ), (31) 



Ai(/ x ,/ 2 )= / dsi / (fe 2 A A (fi,f 2 ), (32) 



is the transverse two-nucleon contraction, ft can be 
seen that the inclusion of NN correlations in nuclei 
leads to a modification of the nuclear thickness function 
Ta^ 3 ) T A l {b). Due to its general structure and the 
basic property J d 3 ri ;2 A^fi, r 2 ) = 0, the sign of the 
contraction is mostly negative, with a small positive con- 
tribution at large separations. In Fig.[T]we present T\(b) 
and AT^(b) for 12 C and 208 Pb. We see that AT^{b) 
is indeed mostly negative, so according to the definition 
Eq. (|3"0|) correlations increase the nuclear thickness func- 
tion and make nuclear medium more opaque [l2j . At the 
same time, the corrections are small, AT^(b) <C T^(b), 
and the effects from higher order correlations, estimated 
in Ref. P3, can safely be disregarded. 

A short description of the way in which the one and 
two-body densities and contractions have been calcu- 
lated is now in order. Following Ref. [l2| , the two- 
body density has been obtained from the fully-correlated 



wave function of Ref. [TH, [3, Vo = F(t)o, where F = 
]li<jEfc=i fk(nj)6 k (ij)} is a correlation operator gen- 
erated by the realistic Argonne V8' interaction [l7|, and 
cf>o a mean field shell model wave function composed of 
Woods-Saxon single particle orbitals. The above wave 
function largely differs from the Jastrow one, featuring 
only central correlations, since the operator F generates 
central {0\ — 1), spin (0 2 (ij) = (fi-ffj), isospin {0${ij) = 
Ti ■ fj), spin- isospin (O^ij) = ■ dj) (f i ■ fj)), tensor 
{Os{ij) = Sij), tensor-isospin (0 6 (ij) = S io (f ■ fj)), 
etc. correlations. The two-body density and contrac- 
tion therefore reflect not only the short range repulsion 
but also the spin-isospin dependence of the interaction, 
particularly that generated by the tensor force. The pa- 
rameters of both the single particle wave functions and 
the various correlation functions have been fixed from the 
ground- state energy calculation so that no free parame- 
ters are present in our approach. 

The contraction A(ri,F 2 ) resulting from our calcula- 
tion exactly satisfies the sum rule J d 3 ri A(r~i, f 2 ) = 0, 
since the one-body density p\ (r\ ) exactly results from the 
integration of p 2 {r\, f 2 ). Notice, moreover, that our one- 
body point density and radii are in agreement with elec- 
tron scattering data [HI]. We have also investigated the 
validity of the approximation in which the nuclear matter 
two-body density p 2 (fi,f 2 ) = px{fi) pi{r 2 ) g{\f\ -r 2 \) is 
used for finite nuclei, finding that it leads to a strong vi- 
olation of the sequential relation J d 3 rp 2 (fi,f 2 ) = pi(fi) 
for nuclei with A < 208. Thus, when such an approx- 
imation is used to introduce correlations in light and 
medium-weight nuclei, a mismatch between the one-body 
density (usually taken from the experimental data) and 
the two-body density is generated. 

Using the nuclear thickness function which includes the 
effects of correlations, Eq. (|3"0"|) , the total cross section, 



Eq. 



acquires a correction, a 



pA 



pA 
°tof 



Aaf£, which 



5 



is positive and can be approximated as, 



Act 



pN 



J d 2 bAT A l (b)exp 



(33) 



which is also positive, since AT^(b) is itself negative. 

We see that this correction to the total cross section 
comes mainly from peripheral collisions, and rises with 
A rather slowly, as A 1 / 3 . Notice that the accuracy of the 
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FIG. 1: Nuclear thickness function T^(b) and the correction 
due to AW-correlations, AT^(b), calculated at the energy of 
HERA-B, for carbon and lead, respectively. 

optical (exponential) approximation in ([8]) is quite good, 
~ 10~ 3 for heavy nuclei, but it gets worse with decreasing 
A, therefore for numerical calculations, as was already 
mentioned, we rely upon the exact Glauber expressions 
throughout the paper. In what follows we neglect the 
real part of the elastic amplitude which gives quite a 
small correction, ~ p 2 p /A 2 / z , and which otherwise can 
be easily implemented. 

The simplest process with a Large Rapidity Gap 
(LRG) is elastic scattering. It is worth noting, how- 
ever, that this channel is enhanced by absorptive correc- 
tions, while other LRG processes considered below are 
suppressed by these corrections. 

The elastic cross section according to (O reads, 



pA _ 



d 2 b 



1 — exp 



1 

'2 a 



M(6) 



(34) 



where T\{b) is given by (pJO")) . 

The quasielastic cross section also gets modifications 
compared to the Glauber expression Eq. (fTB")) . The nu- 
cleon correlations show up in the second order in nuclear 
density, leading to an additional term proportional to 
( a in) 2 AT^(6). Thus, the cross section of quasielastic 
proton- nucleus scattering, pA — * pA* , gets the form, 



pA 

7 qel — 



d 2 b < exp 



-a^T h A {b)-^-AT'X(b) 



exp 



-CTfj(T>) + AT>)) 



(35) 



Notice that in deriving this expression we implicitly used 
the assumption that the impact parameter dependence of 
powers of the amplitude T pN (s) does not depend on the 
power. Although this is certainly not correct, the approx- 
imation is rather accurate as far as the NN interaction 
radius is much smaller than the size of the nucleus. Nev- 
ertheless, we used this approximation only for the sake of 
clarity and simplicity. For numerical calculations, we use 
the more complicated but exact analogue of Eq. 



III. GRIBOV CORRECTIONS VIA 
LIGHT-CONE DIPOLES 

The dipole representation for the amplitude of 
hadronic interactions allows to sum up the Gribov in- 
elastic corrections to all order. We assume the collision 
energy to be high enough to keep the dipole size "frozen" 
by Lorentz time delation during propagation through the 
nucleus. In this limit the calculations are much simpli- 
fied. 

The key ingredients of the approach are the universal 
dipole-nucleon cross section and the light-cone wave func- 
tion of the projectile hadron [5j. Several different models 
were tested in |1G |. by comparing with data on proton 
diffraction. Here we select two models which describe 
diffraction quite well. Both employ the saturated shape 
of the dipole cross section and differ only by modeling 
the proton wave function. 

In the limit of soft interactions the Bjorken x is not 
a proper variable any more, and the dipole cross section 
should depend on energy. We rely on the model proposed 
in [33| and fitted to data, 



crqq(r Tl s) = ct (s) 



1 — exp — 



R 2 (s) 



(36) 



where R {s) = 0.88 fm {s Q /s)°- 14 and s Q = 1000 GeV 2 
[33| . The energy dependent factor ctq(s) is defined as, 



<7o(s) 



'M i + 



3_g|(f) 



(37) 



where (r 2 h )-rr = 0.44 ±0.0 1 fm 2 [l9( is the mean square of 
the pion charge radius. This dipole cross section is nor- 



G 



malized to reproduce the pion-proton total cross section, 



For the proton wave function we employ two models. 



A. q — 2q model 

There are many evidences (although neither of them 
looks decisive) for a strong paring of the u and d valence 
quarks into a small size scalar-isoscalar diquark [2(J l2ll . 
|22| . Neglecting the diquark radius we arrive at a meson- 
type color dipole structure of the proton, 



I* 



N{n,r 2 ,r 3 ) 



exp 



2r| 

p 



(38) 



where ?i are the interquark transverse distances, F3 = 0, 
?t = f\ = fa, and R p is related to the mean charge 
radius squared of the proton as R% = ^§-(r^ h )p. The 
dipole wave function squared, Eq. (|38p . convoluted with 
the dipole cross section, Eq. ([36|) , gives the proton-proton 
total cross section. 

In this model the effect of the Gribov corrections in 
all orders is equivalent to the replacement of the Glauber 
formula Eq. flHJ by, 



i 

altt = 2jd 2 bJdaJ d 2 r T \V N (r T ,a) 



1 _ e -i<rM(rT,B)T%*(b,rT,a) 

-^ g (r T)S )Tl"(b,r T , a ) 



d 2 b 1 



(39) 



Here we consider a qq (or qq — q) dipole of transverse 
separation fx and fractional light-cone momenta a and 
1 — a of the constituents. The integration over these 
variables weighted by the hadron wave function squared 
is denoted as averaging. The new notation Tj[ (6, rr, a) 
is, 



Ti q (b,r T7 



J d 2 l Re T qqN (l,f T ,a)T A (b 



I). 



' a qq (r T ) 

.(40) 

The partial dipole-nucleon elastic amplitude 
ReT qqN (l, rx, a), corresponding to the dipole cross 
section (|36p , was derived recently in [23|, |U, [25| , 



ReT qqN (l,f T ,a) 



0o (s) 



exp 



2 exp 



8nB(s) 
[l + f T (l-a)} 2 ~ 



2B(s) 

[1+ (1/2 - a)f T ] 





(T-f T a) 2 ~ 


+ exp 


2B(s) 



R 2 (s) 



2B{s) 



(41) 



where B(s) = B P J* (s) - |(r^) p - |i^(s). It is easy 
to check that this partial amplitude correctly reproduces 



the dipole-nucleon cross section Eq. (|3"o]) . 

a qq (r T , s) = 2 / d 2 l Rc r** ft r T , a), (42) 

and the slope of the differential elastic pN scattering, 

Bf(°) = 4n l 2 (ReT^ N (l,f T ,a)) . (43) 

These properties of the partial amplitude lead to the 
following relations with the analogous functions defined 
above within the Glauber model, 



vN _ 



'a qq {r T )); 

ReT pN {l) = /Rer^(Z> T ,a<)); 
1 



T'X(b) = 



<*tot X 



<7qq(rT)Tl q (b,r T ,a) 



(44) 
(45) 
(46) 



Thus, the difference between the Glauber formula 
Eq. ijHJ), and the exact expression, Eq. (|3"9")l . is in how 
the averaging over rx and a is done: in the former case 
the averaging is done up in the exponent, while in the 
latter case the whole exponential is averaged. 

Notice that Tj^(b, rx,ct) in the exponent in Eq. (f3"9")) 
can be replaced by T%(b) with a high precision. Indeed, 
at small r\ <C Rq(s) the partial amplitude Eq. (j4"Tj) van- 
ishes as ReT qqN (l, fx, a) oc r T . This rx dependence 
cancels in (|40p with the same behavior of o qq {rx) m the 
denominator. Thus, T q A q {b^ rx, ot) is independent of rx in 
this limit. In the opposite limit of large r\ 3> R\{s) 
the last term in (|39[) vanishes, and the amplitude in- 
tegrated over d 2 l becomes a constant. Moreover, the 
denominator of Eq. (|40[) is independent of rx in this 
limit^ Thus, one can neglect the slow rx dependence 
of T^{b,TT,(x) in Eq. (|46p in comparison with the fast 
varying function a qq {rx), which is equivalent to the re- 
placement T^(b,rx,a) => T^(b). We rely on this ap- 
proximation in Eq. (|39[) and in what follows. 

Thus, for the totalp — A cross section we recover the 
standard expression [a, [13] , 



pA 



= 2 / d l b \ d z r T \^N{rx)Y 



1 — exp 



a qq (rx,s)TH(b) 



(47) 



B. 3g model 

Another extreme is to assume no pairing forces and a 
symmetric valence quark wave function, 



|*Ar(ri, fa, fa)]" 



(Trie)a 



<5(ri + r 2 + r 3 ) 



x exp 



rf + rf + rg 

i?2 



(48) 



7 



Here the mean interquark separation squared is (r 2 ) = 
= 2(r 2 h ) p . In this case one needs a cross section for a 
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FIG. 2: The integral of Eq. (I55(l for carbon and lead, calcu- 
lated at HERA B (dotted curves), RHIC (solid curves) and 
LHC (dashed curves) energies. 

three-quark dipole, 03 9 (ri, r%, f 3 ), where r*j are the trans- 
verse quark separation, with the condition fi+ra+fa = 0. 
In order to avoid the introduction of a new unknown phe- 
nomenological quantity, we express the three-body dipole 
cross section via the conventional dipole cross section <jq q 



CTqq(ri) + <Jqq{r 2 ) + CTqq(r 3 ) . (49) 



This form satisfies the limiting conditions, namely, turns 
into cFq q (r) if one of three separations is zero. 

In this model the Gribov corrections modify the 
Glauber expression Eq. (jSJ) as, 



7 V t i = 2jd 2 bj d 2 n d 2 r 2 d 2 r 3 \V N (n)\ 2 



x < 1 — exp 



-^q{n,s)n(b) 



(50) 



IV. GRIBOV CORRECTIONS TO THE EFFECT 
OF NN CORRELATIONS 



Nucleon correlations lead to further modifications of 
the exponent in Eq. (|39|) . which correspond to the re- 
placement Tl q (b, r T ,a) T q A q (b, r T , a) + AT\ q (b), where 

ATl q (b,r T ,a) = ^- / d 2 l 1 d 2 l 2 Aj i (l ll l 2 ) 

PqqVT) J 

x Rer qq (b~l u r T ,a)Rer qq (b-l 2 ,r T ,a), (51) 



Changing the integration variables d 2 l\d 2 l 2 
where 

L = (h + f 2 )/2; 
5 = h — l 2 ; 



d 2 Ld 2 S, 



(52) 



one has, correspondingly, Aj^(h,l 2 ) A\(L,5). This 
function is rather smooth and varies over distances much 
longer than the interaction radius. Therefore, we can 
take it out of the integral in Eq. (p)Tj) . fixing it at L = 
b. Then, using the partial amplitude Eq. (|41[) one can 
perform the integration over d 2 L in (|5ip and then average 
over tt and a. The result is, 



(53) 



d z L{ Rer qq (h,r T ,a) ReT qq (l 2 ,r T ,a) 

5 2 



pN 



pN 



exp 



AB(s) + R 2 (s)/2 



Here cr^ d is the single diffraction cross section, pN 
XN- and the relation @, 0, 

ReT qq (L,r T ,a 



d l L 



pN 



pN 



(54) 



has been used. 

Data show that at high energies this cross section is 
nearly constant and is about a p s ^ rs 4 mb, the value which 
we use in what follows. 

Notice that data for single diffraction also include the 
contribution from the triple-Pomeron term, which cor- 
responds to diffractive gluon radiation. This term has 
not been included so far in our calculations, which corre- 
spond only to diffractive excitation of the valence quark 
skeleton of the proton (see in [I3|). However, the higher 
Fock components of the light-cone wave function of the 
proton should be also added, which effectively incorpo- 
rate this contribution by using the total single diffraction 
cross section. 

Eq. (fBU)) turns into the Glauber model relation Eq. (f3"Tj) 
if the diffraction term is removed and the denominator 
of the exponent is replaced by ^B P J^ . 

Eventually, the correction related to the nucleon corre- 
lations to the nuclear thickness function, convoluted with 
the dipole cross section, takes the form, 



lA(b) 

I d 2 S exp 



(Tqq(r T ) ATl q (b,r T ,a) 
S 2 



pN 



pN 



AB(s) + R 2 (s)/2 



Ai(S,b). (55) 



The quantity Ia(V) is shown in Fig. [21 for both 12 C and 

208p b 

Since AT a is small and its higher orders are negligible, 
there is no difference between averaging of the exponen- 
tial and of its exponent. Therefore, the result Eq. [55]) 
accounts for all inelastic shadowing effects in A^^V corre- 
lations. Then, the total cross section reads, 



pA _ 



= 2 / d z b 



1 _ e }lA{b) / e -i<74 it ,T^(6)^ 



(56) 
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Here we use a notation which unifies the two models un- 
der consideration. a c ii V is the dipole cross section, and 
averaging corresponds to integration over the light-cone 
momenta of the quarks, weighted with the proton wave 
function squared. 

The results of the calculation of the total, elastic and 
quasi-elastic cross section for several nuclei and HERA 
B, RHIC and LHC energies, obtained within the Glauber 
approach including NN correlations, are presented in Ta- 
bles I, II, III and IV. 



TABLE I: HERA B. 



12 C 


Glauber 


Glauber 


q-2q model 


3q model 






+SRC 


+SRC 


+SRC 


<?tot 


353.71 


364.11 


344.16 


349.37 


Cel 


86.90 


92.96 


82.39 


85.42 


Csd 






5.43 


2.40 


&sd+g 






0.07 


0.06 


O qe 


22.85 


19.62 


21.12 


22.05 


Cqsd 






1.94 


0.84 


<?tsd 






12.47 


12.92 


Cdd 






0.61 


0.26 


27 Al 




Gtot 


697.32 


714.35 


675.93 


688.08 


&el 


201.02 


212.26 


188.22 


196.28 


Csd 






10.82 


4.56 


Csd+g 






0.12 


0.11 


Cqe 


36.39 


31.75 


33.24 


34.86 


Cqsd 






2.92 


1.23 


<?tsd 






19.47 


20.42 


Cdd 






0.91 


0.38 


48 Ti 




<?tot 


1113.52 


1135.53 


1074.67 


1095.93 


Cel 


353.89 


369.77 


327.78 


342.93 


Csd 






16.86 


6.86 


&sd+g 






0.17 


0.16 


Oqe 


48.73 


42.82 


44.57 


46.86 


Cqsd 






3.85 


1.61 


<?tsd 






26.11 


27.45 


Cdd 






1.20 


0.50 



In our calculations nuclear densities which give the cor- 
rect nuclear rms radius have adopted, and this is a rea- 
son of some differences with the results of Ref . [Hj in the 
case of Glauber calculations. The parameters for the 
total nucleon-nucleon cross section and the slope of the 
Glauber profiles have been obtained as in Ref. [10| 



V. DIFFRACTIVE EXCITATION OF THE 
PROTON IN pA COLLISIONS 

While the Glauber model, which is a single channel 
approximation, cannot go beyond elastic scattering, the 
dipole approach treats diagonal and off-diagonal diffrac- 
tive channels on the same footing. Although the calcula- 
tion of exclusive channels of diffractive excitation needs 
knowledge of the light-cone wave function of the final 
state (e.g. see [13, HH), the total cross section of diffrac- 
tive excitation summed over final states is easier to ob- 
tain, since one can employ completeness. Following the 



TABLE II: HERA B (continuation of Table I) 



184 rrr 

w 


Glauber 


Glauber 


q-2q model 


3q model 






+SRC 


+SRC 


+SRC 


(Jtot 


2972.02 


2986.46 


2688.09 


2747.66 


& el 


1174.09 


1187.18 


1025.16 


1074.48 


O sd 






32.27 


11.75 


& sd+g 






0.28 


0.23 


(5 qe 


67.04 


58.35 


57.60 


60.89 


& qsd 






4.92 


2.04 


O tsd 






33.74 


35.66 


0~dd 






1.54 


0.64 


197 An 




(Jtot 


2976.26 


2989.94 


2859.84 


2920.75 


(Jel 


1193.54 


1206.10 


1100.54 


1150.99 


<y s d 






32.55 


11.95 


Usd+g 






0.29 


0.24 


Cqe 


62.94 


54.69 


61.15 


64.53 


0~qsd 






5.31 


2.22 


Ctsd 






35.82 


37.80 


Odd 






1.66 


0.69 


208 Pb 




(Jtot 


3052.11 


3117.62 


2955.57 


3018.21 




1243.00 


1274.60 


1147.01 


1199.14 


Osd 






32.88 


12.02 


&sd+g 






0.29 


0.24 


Uqe 


62.55 


54.11 


61.01 


64.39 


<Jqsd 






5.31 


2.21 


&tsd 






35.73 


37.71 


Cdd 






1.66 


0.69 



standard classification of diffractive channels in terms of 
the triple Regge approach [2§|, one can consider diffrac- 
tive excitation of the valence quark system, which corre- 
sponds to the Pomeron-Pomeron-Reggeon (IP 1PM) term, 
and diffractive gluon radiation corresponding to the triple 
Pomeron term (IP MM). The former mostly contributes 
to small mass excitations, dajdM\ oc 1/My, while the 
latter is responsible for the large mass tail, da/dM^ oc 



l/M x , where Mx is the invariant mass of the produced 
system, pp — > Xp. 



A. Coherent excitation of the projectile valence 
quark system 

The cross section of coherent single diffraction on a nu- 
cleus, caused by excitation of the valence quark skeleton 
without gluon radiation, is given as usual by the disper- 
sion of the distribution of eigen elastic amplitudes, where 
the eigenstates are the dipoles 0, [26[ . 



a sd (pA -> XA) P1PIR = J 



e 2 



(57) 



where I a (b) is given by Eq. {55J) . Dependent on the 
model, the dipole cross section here has the form of ei- 
ther Eq. (pjfj)) , or , and the averaging is weighed by the 
wave function squared having the form of either Eq. (|38p . 
or (1381). 



TABLE III: RHIC. 



12 C 


Glauber 


Glauber 
+SRC 


q-2q model 
+SRC 


3q model 
+SRC 


Otot 


413.71 


425.73 


406.90 


410.20 


<?el 


112.13 


119.68 


109.16 


111.29 


O s d 






3.13 


1.20 


Osd+g 






0.31 


0.30 


Oqe 


26.40 


23.09 


26.13 


26.72 


Oqsd 






0.95 


0.29 


Ptsd 






10.90 


11.14 


Odd 






0.95 


0.29 


208 Pb 




Otot 


3297.56 


3337.57 


3228.11 


3262.58 


O e l 


1368.36 


1398.08 


1314.04 


1343.76 


O s d 






16.78 


5.03 


Osd+g 






1.06 


0.98 


O qe 


66.06 


58.47 


71.99 


73.92 


Oqsd 






2.39 


0.56 


Otsd 






30.03 


30.83 


Odd 






2.39 


0.56 



Although Gribov corrections to the total cross section 
are known to be small, well within 10% [3(| HH, this is 
because they affect only the second exponential term in 
Eq. ([5]), which is small. However, this term itself is mod- 
ified significantly by the inelastic shadowing corrections. 
Therefore, one should expect a considerable increase of 
both terms in (|57|) due to inelastic corrections, which 
make the nuclear medium considerably more transparent 



compared to the Glauber model Nevertheless, both 
terms are small for heavy nuclei and suppress diffraction 
everywhere except at the nuclear periphery. Thus, the 
cross section of single diffraction should rise as A 1 / 3 , with 
a coefficient which is sensitive to the inelastic shadowing 
corrections and NN correlations. 

The details of the calculations with both models under 
consideration can be found in [lfj . The numerical results 
for several nuclei and energies are presented in Tables I, 
II, III and IV. 

TABLE IV: LHC. 



12 C 


Glauber 


Glauber 


q-2q model 


3q model 






+SRC 


+SRC 


+SRC 


Otot 


598.79 


613.68 


591.05 


592.12 


Oel 


198.11 


208-59 


194.84 


195.65 


Osd 






0.74 


0.20 


Osd+g 






2.58 


2.56 


Oqe 


49.10 


45.42 


45.03 


45.22 


Oqsd 






-0.66 


-0.86 


Otsd 






6.97 


7.00 


Odd 






-0.66 


-0.86 


208 Pb 




Otot 


3850.63 


3885.77 


3833.26 


3839.26 


Oel 


1664.76 


1690.48 


1655.70 


1660.67 


Osd 






2.62 


0.59 


Osd+g 






2.58 


2.56 


Oqe 


120.92 


112.65 


113.37 


113.88 


Oqsd 






-2.08 


-2.62 


Otsd 






17.55 


17.63 


Odd 






-2.08 


-2.62 



B. Coherent diffractive gluon radiation 

Diffractive gluon radiation also contributes to the sin- 
gle diffractive process pA — » XA. Correspondingly, the 
single-diffraction cross section Eq. (f57|) must be corrected 
for this excitation channel. The cross section of coherent 
gluon radiation on a nucleus is given by [33| . 



<j sd (pA XA) 3P = — In 

4-7T 



s(l - Xq) 



d 2 be u W 



1 — exp 



e 2 



1 _ rph 

2 °dip J A 



Ml 



d 2 r T 



9 

16 



l {rr,8)T^b)-%I A {),) 



58) 



Here Mq = 5 GeV 2 is the minimal effective mass squared 
of the proton excitation, xo = 0.85 is the minimal value of 
Feynman x, which can be treated as being in the domain 
of the triple- Regge kinematics. [29| . 
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The first factor in (|58|) accounts for the absorptive cor- 
rections, which are due to the lack of initial/final state 
interaction of the valence quarks propagating through 
the nucleus. Further details about the calculations can 
be found in [101 ] . The numerical results for several nuclei 
and energies are presented in Tables I, II, III and IV. 



VI. QUASIELASTIC SCATTERING WITH AND 
WITHOUT EXCITATION OF THE PROJECTILE 

In the cases when either the beam proton (pA — > XA) , 
or the nucleus (pA — > pA*), or both (pA — > XA*), are 
diffractively excited, one can make use of completeness, 
which substantially simplifies the calculations. As was al- 
ready mentioned, the important condition for the nucleus 
is that it decays into nuclear fragments with no new par- 
ticle produced. The dipole formalism for these processes 
was developed in [l(| ■ Here we rely on those results and 
introduce corrections related to NN correlations. 

The simplest processes are double excitation pA — > 
XA*, where X includes the ground state proton, as 
well as quasi-diffraction, with all diffractive excitation of 
its valence quark system (without gluon radiation), and 
breakup of the nucleus. All channels of coherent inter- 
actions which leave the nucleus intact, should be sub- 
tracted, 



separately for both the incoming and outgoing protons, 



VqeiipA -» pA*) + a qsd (pA 



d 2 b I e~°^™ 



■+ XA*) 

>li v T A^ 



k=0 



1 


\ n(b) - 


k 


e /x(6) 











jA(b) 



>k0 



)2k 



d{T'X) 2k 



-o dlv T H A {b) 



(59) 



Here besides the function 7a (fr) defined in (|55p we intro- 
duce a new one, 



h(b) = j d 2 hd%A\(h,i 2 ) 

2 Re r« a (fi, f T , a) - (Re T^{h,f T , a) 
2 Re (h,f T ,a)-( Re r« 9 (f 2 ,r T ,t 



pN pN 



'tot 



pN 
sd 



-1 2 



vN 



lA(b) 



(60) 



In order to simplify the calculations, we neglect here the 
difference in the slopes of powers of the partial amplitude. 
This is a second order correction, i.e. a correction to a 
correction. 

One can single out in (159")) the quasielastic channel. For 
that purpose one should average over the dipole sizes, 



pA _ 



'qel 



= d'b 



e lA(b) e - 



,n(b) 



jA(b) 



d 2 b 



E 

k=Q 



1 





k 







jA{b) _ e lA(b) S 



ki ) 



(61) 



This is a fast converging series due to the smallness of 
the elastic cross section. We control the accuracy to be 
within 1%. 

Subtracting (|6T|) from ([59]) one can get the quasi- 
diffractive cross section, which includes the proton ex- 
citations without gluon radiation. To include gluon radi- 
ation we use the same prescription as in (|54[) . replacing 
the TP IPM term, [o'^] ]P1 p 1R , by the total single-diffraction 
cross section. 

In the case of nuclear breakup the recoil bound nu- 
cleon can be also diffractively excited. We relate the 
cross sections for such channels to the above calculated 
quasi-elastic and quasi-diffractive processes, in the same 
way as in [ic| . 



VII. GLUON SHADOWING 

In terms of the parton model, gluon shadowing is inter- 
preted in the nuclear infinite momentum frame as a result 
of fusion of gluons originating from different bound nu- 
cleons. This process leads to a reduction of the gluon 
density in the nucleus at small x. The ultimate form of 
gluon shadowing is gluon saturation [32] . 

In terms of the dipole approach gluon shadowing is 
described as Glauber shadowing for higher Fock states 
containing gluons, [33| . The effect turns out to be rather 
weak due to the shortness of the quark-gluon and gluon- 
gluon correlation radius, an observation which is sup- 
ported by many experimental evidences [3J, |35j . For this 
reason we neglect the small effects of nucleon correlations 
in the calculation of gluon shadowing, and use the results 
of Ref. 0. 



VIII. CONCLUSIONS 

In this pap er we further developed the dipole approach 
of [H, llOl ] to the calculation of Gribov inelastic cor- 
rections. We employed two models for the proton wave 
function, which result in reasonable diffractive cross sec- 
tions for pp collisions. Here we increased the accuracy of 
the calculation of the cross sections of different diffrac- 
tive processes on nuclei by improving the model for the 
nuclear wave function. Namely, we went beyond the pop- 
ular single particle density approximation and introduced 
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corrections for nucleon-nucleon correlations, which lead 
to sizable effects, modifying the effective nuclear thick- 
ness function [l2j • While inelastic shadowing corrections 
make the nuclear medium more transparent for colorless 
hadrons, the nucleon short range correlations work in the 
opposite direction making the medium more opaque. The 
influence of both effects on different diffractive channels 
vary. They are especially large for quasielastic and quasi- 
single diffractive processes associated with the survival 
probability of colorless hadrons propagating through a 
nuclear medium. 
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